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Abstract 

Itoyama-Tokura type USp matrix model is discussed. Non-Abelian Berry's phases in a T-dualized model of IT model were 
reconsidered. These phases describe the higher dimensional monopoles; Yang monopole and nine-dimensional monopole. They 
are described by the connections of the BPST instanton on S 4 and the Tchrakian-GKS instanton on S 8 , respectively. 

As a preparation to understand their effect in original zero-dimensional model, we consider partition function of simplified 
matrix models. We compute partition functions of SU, SO and USp reduced matrix models. Groups SO and USp appear in 
low energy effective theories of string against orientifold background. In this evaluation we chose different poles from that of 
i Moore-Nekrasov-Shatashvili and our previous result. 

. The position of poles explain branes' and the orientifold's configurations. There is a brane which is sitting on the orientifold 
in the SO(2iV) model, while in USp(2A r ) and SO(27V + I) model there are no branes on the orientifold. The grand partition 
functions of these models are considered. They follow to linear second order ordinary differential equations and their singularities 
are q = 0, do. Their solutions can be analytically continued to whole q plane. We show the expectation values of the number TV 
of A and C cases as examples. There is an ambiguity coming from the problem on sign. Grand partition functions with minus 
sign give effective actions which have cusp singularities. 
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I. INTRODUCTION 



In 1997, Itoyama and Tokura considered a matrix model which we call USp matrix model. The model was studied 
in order to understand the dynamics of string against the orientifold background pj . There is another USp matrix 
model by S.J. Rey and N.G. Kim Q and their work must be important for us. Because the eleven-dimensional 
minimal coupling supcrgravity theory contains rank-three antisymmetric tensor fields, the theory has been considered 
as the low-energy effective theory of the extended object "memferane" whose trajectories of movement are, in general, 
2 + 1-dimensional space. The mass spectrum of the supcrsymmctric membrane is described by the Hamiltonian of 
the dimensionally reduced model of the ten-dimensional supersymmetric Yang-Mills theory to one dimension. 

Type IIA superstring theory has infinitely many BPS states which couple with the Ramond-Ramond one-form and 
these BPS states are identified with configurations of D-particles. The tower consisting of these configurations of 
D-particles can be interpreted into the Kaluza-Klein modes accompanied by the circular compactification along the 
eleventh direction. This viewpoint leads people to the discussion on the duality between the theory of membrane and 
IIA superstring theory [H, Q . The matrix theory has the time direction and it is not clear whether the model has 
covariancc or not. Ishibashi, Kawai, Kitazawa, and Tsuchiya have proposed a matrix model whose masslcss spectrum 
is the same as that of type IIB string theory It was conjectured that the large- TV reduced model of ten-dimensional 
super Yang-Mills theory can be regarded as a constructive definition of string theory. 

Z=jr f [dXWd*W] exp(-SW) , S N = a [-^[Jf^, Jfflf - Wr*^, &>0]\ + /?Trl , 

where xf ] and tf W are bosonic and fermionic N x N Hermitian matrices, respectively and p, = 0, 1, • ■ ■ ,9. a = l / g 2 
is the inverse of string tension or gauge coupling and f3 is the chemical potential which is needed for insertion of one 
instanton. 

The IKKT matrix model was constructed as the united theory of string. If in fact the model unifies string theories, 
the model should not depend on the background. The background independence is the most significant implication in 
the inclusion of quantum gravity in string theory. The geometry of space-time should not be set up a priori, rather it is 
generated by a highly nonperturbativc effect, the condensation of strings. Therefore understanding of the background 
independence is promisingly the key ingredient to seek the underlying principle of nonperturbative string theory. 

We can include another type of background. Oricntifolds are generalized orbifolds. In the orbifold construction, 
discrete internal symmetries of the world-sheet theory are gauged. In the orientifold, products of internal symmetries 
with world-sheet parity reversal are also gauged. Roughly speaking, these symmetries yield restrictions of gauge 
groups to unitary symplectic (USp) or orthogonal (SO) groups. 

In this article, we will consider partition functions of reduced matrix models. Several people evaluate the quantity. 
The quantity relates to the Witten index [|| of matrix quantum mechanics and to the dynamics of D-particles in the 
context of M-theory H, H, foL flQ| ■ Monte-Carlo simulation has been considered [ill, Ell [HI and the convergence of the 
integral was studied in [lj] . Moore-Nekrasov-Shatashvili obtained the result 1 /N 2 by using the deformation method 
and the Monte-Carlo result by Krauth-Staudacher [ll[ agree with the MNS result. We start from the MNS's contour 
integral and obtain the different result from the result 1/N 2 for SU(N) case. There are four types of the partition 
functions Zrmm, ^mc, Zmns, Zx- The true partition function of the reduced matrix model with suitable normalization 
factor is denoted Zrmm, the result of the Monte-Carlo simulation which might be used different normalization factor 
is Zmc, the MNS result is Zmns and ours is Zx- The partition function Zmc is a finite list, however it agrees with 
■Zmns f° r SU(2) and SU(3) cases. Let us start from doubt of the equality between Z-^mm and Zmns- It implies that 
we do not believe the equivalence between Zrmm and Zmc- For SU(2) case, our result agrees with their result 1/4, 
while SU(3) case there exist the difference of factor 2. Two quantities -Zmns, Zx are different with each other. The 
difference is occurred from the difference of the choice of the contour. In this paper, we will not discuss the justice 
of these results, but will discuss the meaning of the multiplicity. Our resultant effective potential for SU series has 
slightly similar shape with a graph in [l5j . however our independent variable is the expectation value of the matrix 
size and their variable is temperature. 

There are two kinds of bosonic sector of Type IIB superstring theory; i) NS-NS sector: dilaton <fi, graviton g^o, 
two- form field Bj,~- ii) R-R sector: axion y, two- form field self-dual four-form field Af--~. The combination of 

these scalar fields r = \ + i cx p(~</ ) ) 1S regarded as a modulus of torus fiber T 2 . Then the doublet (i?J 5 , -B£>) becomes 
a SX(2,Z) multiplet. Therefore IIB superstring on a compactified space B is considered as a compactification of F 
theory over M, where fiber bundle ir : M — > B (tt^ 1 (z) ~ T 2 ) is defined by the moduli function t(z). One famous 
example of this configuration is M = A'3 and B = CP 1 ~ S 2 . In this case, F-theory over A3 is dual to IIB over S 2 . 
Here K3 is smooth Calabi-Yau 2-fold. Especially F theory over A3 near the orbifold limit of A"3 describes the dynamics 
of orientifold of Type IIB over T 2 . On the other hand, orientifold of Type IIB on T 2 is T-dual to Type I on T 2 . One 
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Realization of A'3 (not generic) is a quartic surface in CP 3 ; V(f) = {Z G CP 3 |/(Z) = 0, deg/ = 4, / : homogeneous }. 
The cohomology class is obtained by the Kahler 2-form on CP 3 ; (1 + J) 4 / (1 + 4 J) ~ 1 + 6 J 2 . From this, we can read 
off that the Euler number is equal to 24. 

Let us talk about Itoyama-Tokura Model which is a kind of reduced matrix model. The model is written in terms 
of dimensional reduction of the d = 4, Af = 2 USp(2A r ) supersymmetric gauge theory to zero dimension. The model 
has i) one adjoint vectormultiplet; ii) one anti-symmetric hypcrmultiplct; iii) Nf fundamental hypermultiplets. The 
total action of this model consists of two parts; Srot — Sc + So ■ Here the closed sector Sc is obtained by orientifold 
projection from IKKT action and the open sector So is from space-time filling D-brane. Three representations of 
USp group appear. The adjoint and the antisymmetric representations are, in a sense, parts of Hcrmitian matrices; 

u(27V) = {iX G M(2N,C)\X^ = X} , usp(2N) = {iX G u(2N)\X T J + JX = 0} 

asym(2N) = {iX G u(2N)\X T J - JX = 0} , J = ( ° T T * 

y—iN u 

Therefore u(2N) is a direct sum of adjoint and antisymmetric; u(2N) ~ usp(2N) © asym(2N). Let iX be an element 
of usp(2N) and iY be an element of asym(2N). The block notation show us their construction 

( A B \ (CD 
X - \Bi -A T ) ' Y - [pi C T 

where At = A, C+ = C, B T = B, D T = -D. 

The orientifold projection is described as follows. The Af = 4 SU(2iV) supersymmetric Yang-Mills theory consists 
of four multiplets V, $i,$2,^3 which are all adjoint representation of SU(2./V). Let us project them into Af = 2 
USp(2A r ) supersymmetric gauge theory with adjoint representations V, 3>i and antisymmetric representations <J> 2 , $3. 

We also add Nf fundamental representations Qi,Qi by hand. In other words, a ten-dimensional vector splits into 
six-dimensional vector and four scalars; X^ G su(2N) =>- X^ G usp(2N) (/i = 0, 1, • ■ ■ ,4,7) © X a G asym(2N) 
(a = 5,6,8,9). Splitting of a Majorana-Weyl fermion ^ which belongs to the SU(27V) adjoint representation can 
be considered; $ = ip + A, where ^/i is a USp(2A r ) adjoint eight component spinor and A is an anti-symmetric eight 
component spinor. The action of the closed sector, Sc, is 

So = —TrlX^Xp] 3 - ^Tr*rA[X A) tf]. 

Let us omit the explanation of the open sector So- In this model, the classical configuration is given by diagonal 
matrices Xm and all fermions are set to zero; \& = 0, Q, Q = 0. 

Xn =diag(x^, ••• ,x%,-x1,--- ,-x*) , X a = diag(xi, • • • ,xf , a?i, • • • ,x%) 

They form TV pairs of ten-dimensional vectors; {(a:*, a;*), (— x*, x l a )}i=i,.-- ,N- In each pair, one vector is a mirror 
image of the other with respect to the orientifold. In order to reveal the physics of the fermion of the fundamental 
matter, we computed the non-Abelian Berry's phases against the classical background [3, [I3|- ^ was computed as 
the fcrmionic integration of the matrix model. In the paper, we did not consider the effect from the fermionic diagonal 
terms. In order to make the discussion precise, we should take account into the bilinear terms with respect to the 
fermionic diagonal elements, which are yielded from the quadratic completion of fermionic off-diagonal terms. In this 
paper we will not treat this problem. In fact such a term has influence on the integration of the bosonic off-diagonal 
degrees. Suppose that X^ is diagonal. We split the fermion into diagonal and off-diagonal terms; ^ = + ^o- As 
we mentioned above, we neglect ^d- Fcrmionic part of the whole action is 

S F = -^Trf r^(acLY A )*o + £ {xPAi^ + Mi)xP + X^Jil^ + M^} , 
^ i=l 

where 7 M are six-dimensional Dirac matrices and the eight-component spinors 1 xf 1 " 1 are components of the 
fundamental hypermultiplets Qi,Qi, respectively. The nf x rif matrix M. is the mass matrix of the fundamental 
matter and please permit to restrict our discussion in the case that the mass matrix M is diagonal. The off-diagonal 
part can be written in the components using roots and weights; 
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where A is the root system of usp(2N) and A' is a set of non-zero weights with respect to the anti-symmetric 
representation asym(2N). T a and T w are generators of usp(2N) and weight vector of asym(2N), respectively. They 
are represented as unitary matrices. Let H = diag(/ii, • ■ ■ , h^, —hi, ■ ■ ■ , —Hn) be an element of a Cartan subalgebra. 
Then the dual basis is defined as ei(H) = hi (i = 1, 2, • • • , TV). 

A = {±(ej - e j ),±(e i + e,-), i < j, 2e h i =!,■■■ , TV} 
A' = {±(e 4 - ej),±(ei + < j} 

Components ip a , A w from roots and weights which are explained like ±(ej — e-,), ±(e 4 + e^) form sixteen- 
componcnt spinors and components from ±2e^ form eight component spinors. In other words, the composition 
(V'iCe,— e )j -^±(ei-e )) an d (ip±r ei +e •) i ^±(e»+e ■)) form sixteen-component spinors, while (V^eijO) are still eight- 
component spinors. Therefore we obtain two types of action; i) <S7[£, z] — £7 M Zp£ where £ is a 8 component spinor, 
and z is a 6-dimcnsional vector; ii) Sj/[3, y] = ET^y^S where S is a 16 component spinor and y is a 10-dimcnsional 
vector. However this model is a model on zero-dimensional space (a set of discrete points). Let us consider T-duality 
in order to "make" time direction and let us concentrate on the "one-particle state" . Then Hamiltonian becomes 

ff 7 [z]= 7 %, HnM= T %- 

We set yo = and zq = and evaluate their Berry's phases in consideration of degeneracy. These Hamiltonians 
are higher dimensional generalizations of the system where spin couple to background magnetic field. In the case of 
Hi [z] , the background magnetic field are given by the Hodge dual of four form field strength and the corresponding 
gauge fields are three form field. Because of the degeneracy of these Hamiltonians, we obtain non-Abelian connections 
which are related to the generalized Monopoles. They are generalization of the Dirac monopole in five-dimensional 
space and nine-dimensional space. Berry's connections are 

z = {zi,z 2 , z 3 ,Z4,z 5 y , y = {yi,--- ,z 9 y , R 2 Z = z 1 ■ z , R 2 = y f ■ y , 



1 / T 2 

Z = - {z 4 + i{z\G\ + z 2 a 2 + z 3 a 3 )} , T=\jR% — z\, A Yang = dZZ 1 , A = R z + z 5 . (1) 

T » T- + A" 

Hypcrsurfaccs defined by the condition A = const are hypcrboloid and the hypersurfaces shrink to the Dirac string 
after taking limit A — ► 0. This five-dimensional monopole is called Yang monopole [l8j . If we transform this space 
to a space where the hypersurfaces A = const, become hyperplanes. There appears a mirage four plane where A = 0. 
The metric on the hyperplanes, which is induced from Euclidean metric on the original five-dimensional space, has 
non-trivial curvature. In fact, such a hypcrboloid is a curved space. We omit the exhibition of the Tchrakian-GKS 
type solution. Thus the Hamiltonian Hi gives BPST instanton connection on S and the connection satisfies the 
self-dual equation; F = ± * F. While remaining Hamiltonian Hn gives the Tchrakian or GKS connection on S 8 and 
the connection satisfies the generalized self-dual equation F/\F = ±*FAF [TH, [2(| . They are generalizations of Dirac 
monopole and the latter is generalization of the self duality. In nine dimensional space, the shape of singularities look 
like a point. These brancs which arc denoted by solid lines come from fundamental multiplets and their distances are 
determined by the mass matrix. 



a • • D-particle 



.1 

ijif 



TV/ D4-brane 
Orientifold 4-plane 



FIG. 1: D-particles are from the Hamiltonian Hn and orientifold and D4-branes are from Hi. 



Our notation in [T3| owed to [2l| which was introduced by Itoyama. The non-Abelian generalization of Berry and 
Simon (22l, l23j is argued by Wilczek and Zee [24| . The derivation of BPST instanton as Berry's phase is discussed in 
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II. TCHRAKIAN'S MONOPOLE 



Because the Yang monopole is a singular object like Dirac monopole, we cannot determine its mass or energy. 
In this section we will consider the Tchrakian's monopole 0, [26| as a possibility of the regularization of the Yang 
monopole. He constructed five-dimensional finite energy monopole solution which is an analogue of the 't Hooft- 
Polyakov monopole. The computation of the Berry's connection was done with adiabatic approximation. If the 
parameter z is near the origin, the adiabatic approximation is not valid because the gap of the spectrum of Hj is not 
large enough. Such a gap is needed to avoid the transition between different energy eigen states. Therefore in order 
to clarify the behavior around the singular point, we need such a regularization. 

Let us work on the five-dimensional Euclidean Clifford algebra: {7a, 7b} = 2S a b, where a, b = 1,2,3,4,5 and 
(Tab = [7o)7fe]/2 are SO(5) generators. These gauge fields are brought together in one differential form which takes 
value in the Lie algebra; A = (l/2)A h ^abcdx a . The scalar fields are also represented as a matrix <j> = <p a ^ a - Field 
strength two form is F = dA + gA 2 where g is a gauge coupling constant. The energy is defined as 

E = J Tr jg-^jCF A F) A *(F A F) + ^D<j> A *D<f> + XV (<j))d 5 x 

The symbol * is the Hodge dual operator with respect to the Euclidean metric on R 5 . From this energy, the Bo- 
gomol'nyi equation with the Prasad-Sommerfeld limit becomes D<f) = ± * (F A F). Let us consider the Hedge- Hog 
solution; 

e = x a la , r 2 = x a x a , A = 1 ~ K ^K de, <f> = H a U{r) ei 

where Hq is the vacuum expectation value. Unknown functions are K(r) and U(r) and their boundary conditions are 
U(0) = 0, K(0) = 1, U{oo) = ±1, K(oo) = 0. After change of the variable a 3 = 2g 2 H /3 : s = ln(ar), X(s) = K 2 , the 
equation becomes an autonomous differential equation: 

Y(X) = —r -— , X(l- X)Y— =2XY 2 + 3Y-2 . 

V ' X(l - X) ds ' V ' dX 

This equation is the Abel differential equation of the second kind. Let us put Z = dY/dX. Points (X,Y,Z) 
are sitting on a surface defined by a quartic polynomial. The surface is singular and has a line singularity. The 
differential equation does not have the same property as that of Kovalcvskaja. Fortunately numerical evaluation 
shows the existence of a flow which connect two boundary points. The effective theory against this background are 
written in terms of three form gauge field. 



III. PARTITION FUNCTION OF LITTLE SO AND USP MATRIX MODEL 



In this part we will compute the partition function of matrix models whose actions are given by the dimensional 
reduction of d = 4 Af = 1 supersymmetric Yang-Mills theories. Here gauge groups are classical groups and we list some 
properties of Lie groups in the appendix. Low dimensional counterparts of IKKT matrix model have been discussed 
in [27] and more. We use the terminology "litte" matrix model borrowed from them. The partition function of matrix 
models are considered in the context of the discussion about the existence of D- particle bound states; @, H, H, fiot - 
They evaluated the Witten index of matrix quantum mechanics @. We will follow to the equivariant deformation 
method which was used in Q. Once we performed the computation [28| . In this paper we will obtain different result 
from them. 

The partition function of a reduced matrix model is given as a matrix integral. Moore- Nekrasov-Shatashvili obtained 
the following result for d = 10 Af = 1 supersymmetric Yang-Mills theory. 

v ' k\N 

Here the summation with respect to k is taken over the set which consists of all positive factors of N. In order to 
obtain these result, Moore-Nekrasov-Shatashvili used equivariant deformation method. 

Recently we have applied the deformation method to Af = 1, d — 4 in the case of SO and USp gauge theories. 
In this paper, we retry the calculation and will study their grand partition functions. This work is a preparation to 
evaluate the integral ob the USp matrix model. 



5 



Q <j> (f> 7] x B Vt Tpi 

12-2-1-100 1 
TABLE I: The assignment of the ghost charge. 

Let us consider a model whose action is 

1 r ia 1 



S[As-i] = Tr ^- — [ Vm ,v n Y - ^A7"> m , A] 

Here the gauge group is SU(AT)/Zjv and suppose that the metric is Lorentzian to, n = 0, 1, 2, 3. This action is a low 
dimensional counterpart of the IKKT action. The signature of metric should cause a problem on its convergence. 
Suppose that the integration path of <j> is rotated. We will follow to the MNS's method and the matrix integral 
reduce to the residue calculus. We will skip the detail of the derivation. The degrees of freedom are V m , A which are 
N x N Hcrmitian bosonic and fcrmionic matrices, respectively. Let {T a } a=1 ... i j\r 2 -i be generators of SU(AT). The 
matrices are expanded by these generators: v m — v mTa, A = ^ A°T a . The components A° are Majorana spinors: 
A a = (A a , A a )*. These are not fields, but constant matrices. Let us rearrange these matrices; <\> = —(vq + V3) ,<f> = 
vq — V3 , A 1 = {ipi + iip2)/V% , A 2 = r)/2\/2 + i\. Let us use an auxiliary field B and modify the action S — > S[B\. 

1 fi 2 1 

- ^2 Tr I - J2 ( & vA [0, v r ] - V/ Vj] + V/ [«/ M ) ~ g [0, 0] 2 (2) 

- J B 2 + i £S + ir,[0,7 / ]+ X [0, X ]+V2 X ([ Ul ,^ 2 ] + [^ 1 , U2 ])l . (3) 



Here £ = y/2i[vi,v 2 ]. A BRST-like charge Q is defined as Q<f> = 0, Qcj) = £ , = [0, </>]_, Qv/ = V/ ,2^/ = 
[(/), uj] , Qx = B , QB = [(f), x) j where Q is nilpotent up to traceless part: Q 2 (j> = , Q 2 <fi = [4>, <p], Q 2 vj = [(f), vj] , etc 
(I = 1,2). The action S[B] is written in the Q closed form, 

s[b] = qv, v = j 1 2 <^y,M^] + iv[M] + xB-ix£\ ■ (4) 



kg 

Thus, S[B] is Q-exact. The partition function is not a path integral, but an ordinary integration, 

z = ^J^j J [dv] [dX] [dB] cx p(^ p ) • ( 5 ) 

All fcrmionic contents form a finite-dimensional Grassmann algebra which are generated by {A^}. The Grassmann 
algebra are graded and we call the number of multiplied A° rank. Let us understand that the fermionic integration 
pick up the highest rank elements. This integral might be divergent because there are flat directions, though we take 
account into the Euclidcanization of (f> and there is the vanishing determinant term which comes from the fcrmionic 
part. Let us change our space to Euclidean space —S[B]e '■= iS[B]. 

Ze = J W>] m [dB] [dvrlidr]} [d X ] exp(-QP) = -J-- J [#] cxp(-5 off ) . (6) 

Let us assign ghost charge +1 to Q. The Lagrangian should have ghost charge 0. If the ghost charge of \ 1S 
v, the auxiliary field B has ghost charge u + 1. Then the term QTr^i? should have charge and it means that 
1 + v + (u + 1) = 0. So v = —1. The matrix <f> have the same charge as that of Q 2 . The charge of <f> is 2 and <f> has 
—2. From the term QTrx£ , we can read off the charge of vj as 0. Let us deform the BRST charge with respect to 
the little group SO (2): Q E , 

QeVl =1>i, Qs4> = t , QeX = B , Q s (f> = , 

Q £ Vi = [0, vj] + euEvj , Q e £ = [0, $\ , Q e B = [(f>, X ]+E x . 

Here E is the deformation parameter and let us consider the deformed action: S E = Q e V . Let us assign the ghost 
charge of E as 0. Then the additional terms coming from this deformation have ghost charge —1 and it does not 
matter. Therefore we can consider this deformation does not change the value of integration. 



6 



Localization technique shows that the integral reduces to counting of the fixed points. The integral is Gaussian for 
almost all fermions and B and <j>. We add mass terms for the convergence of zero-modes' integrations by hand. After 
that in the classical limit (this procedure picks up the fixed point), we can integrate them and obtain the simple result 

1 'w, 1 



vol(G) J 1 rj Detadj.(ad<?i + -B) 

where we neglect various numerical factors coming from the Gaussian integrations. We expect that such a numerical 
factor cancels for pairs of boson and fcrmion except for <j>. In addition, though the first integral converges, the resultant 
integral diverges. It means that there exists a space to discuss the validity of the deformation method. Let us consider 
this integral with some regularization, which may not be equal to the original integral. The integrand is invariant 
under a transformation (f> — > <fi' — UtpU" 1 with a unitary matrix U . Let us use the Wcyl integration formula. These 
integrals reduce to integrals over the Cartan subalgebra. The volume of group is given as vol(SU(7V)/Zjv). The 
integral can be reduced to a simple integral over diagonal traceless matrices. 

where N in the numerator is order of the center of SU(iV) and N\ is the order of its Weyl group. Generalizations 
to other groups are obtained immediately. Let $ be one of root systems A^-i, Bn,Cn, Dn- The integral for $ is 
deduced as 

* Z jfr^iTT M) 



Here in order to obtain finite result, we close the contour with additional paths. We will talk the detail later. The 
symbol a denotes a root, the variable <j> is an element of CSA, Z is the center of the universal covering group of G, and 
W is the Weyl group. This integral diverges because poles sitting on the contour. In order to avoid such a divergence, 
let us shift the parameter E to the imaginary direction Im(E) > 0. Then we count poles where (a, fa) > for all 
positive roots a because we chose the upper contour on each complex fa-plane. This region is a Weyl chamber. Let 
us define quantities pi = (sj, fa) where Sj(£ = 1, ■ • • , r) arc simple roots. The linear transformation <fi — > (pi) yields a 
Jacobian 1/G where G is the determinant of the Cartan matrix and G = Our previous work shows that each 

nontrivial pole which exists inside of a Weyl chamber is given by pi — E for all i for all classical gauge groups, where 
Pi are defined by the fundamental root system obtained by the Weyl chamber. Let us put ourselves on a skeptical 
ground. In this article we do not treat the result by MNS as exact. Our selection of contour excludes the additional 
multiplicity. This is not a truncation. This differs from that of MNS. This contour is, in a sense, a fundamental cycle 
and MNS count {N — 1)! cycles. Remaining poles are sitting on boundaries of Weyl chambers. The consideration on 
poles on boundaries are more sensitive. In this paper we ignore such a pole on boundaries. Such a poles might exist in 
SO and USp cases, while in SU model, there are no poles on boundaries. Concrete calculation show that the residue 
at poles on boundaries in some cases in SO and USp models vanish. This is the reason why we do not consider the 
contribution of poles on boundaries 29]. In addition, diagrams which are not connected give points on boundaries. 

Every root a is written as a linear sum of simple roots with integral coefficients and this fact inherits to the variable 
{a, fa). 

r r 

a = ^u(a)iSi , (a, 0) = }] u(a)iPi . (7) 

i=l i=l 

Therefore the integral is written as 



zm = _L rr rr HWj n u ^ 

#W \ Al ll u ( a )._i I \ 11 u(a) l + l 



In order to obtain the result 1/N 2 , we need multiplicity (N — 1)1 for SU(iV) cases. We do not insert the multiplicity 
in this paper. We show several value for low rank groups in Table [TTJ These poles are written as 

An-i ■ 4>% = {N + l)/2 — i, Bjf.<fH = N-i + l, C N : fa = N - i + 1/2, D N : fa = N - i, (8) 
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Ax 


Ci 


Ai ffiAi 


L>2 


B 2 


C 2 


A3 


£>3 


1/4 


1/4 


1/4 x 1/4 


1/16 


3/64 


3/64 


1/96 


1/96 



TABLE II: This table show that the resulting residues correspond with the isomorphic algebra. 



where 1 < i < N. We can interpret these poles as a configuration of branes. Residues arc computed for all classical 
groups. 

1 1 (2N -3)!! 1 (2N-1)\\ 

Z \ A N-i\ = > z Wn) = 2 n n . ni {2N _ 4)! , > z \ b n] = Z[C N ] = 2 N + i N . Nl{2N _ 2)u > ( 9 ) 

Equivalence relations, Z[£?tv] = Z[Cjv] and (N + 1) 2 Z[Dm+\] = NZ[Cn\, are derived. 



A. Grand Partition Functions; A,C 



Now we can sum up them. Let us define corresponding grand partition functions; 



Q[A;q] := 1 + £ Z[A N ^]q N = 1+ £ 

oo oo 

Q[D;q] :=l+J2 Z [ D N]q 2N = l+Yl 



N=2 



N=2 



N M 

1 (27V -3)!! 1 (q 
N (2/V-4)!! AH 



6[B; g] := g + £ Z[ J B w ]g 2JV+1 = g + g £ 



(2JV-1)!! 1 fq 2 



N=l 

oo 



/V 



e[C;«] :=l+^Z[ J B Ar ]g 2W = l+^ 



^ (2JV)!! TV! 

(27V - 1)!! 1 ^g 
(27V)!! TV! 



2V 



AT 



(10) 

(11) 

(12) 
(13) 



where g = e _/3 and /? is the chemical potential which couples to the matrix size TV. The convergent radii of these 
functions are infinite. These functions satisfy the relation 0[£>; q] = gO[C; g] and written in terms of the exponential 
integral function Ei(z) and the confluent hypergeometric function ^>(s,t;z) defined in the appendix. Here we study 
about grand partition functions of A and C series; 



e[A;q]. 
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-dx = Ei(g) C ln(g) - g + 1 



6[C;g] + = $(l/2,l;g 2 /2) , (14) 



These functions, 0[A;g] + and 0[C;g]+, obey the following differential equations. The meaning of the subscript + is 
explained later. 



d 2 _ x d 

df- {1 - q ] dq-. 



&[A;q}- 



1 



d 2 _ x d 

d^' iq - q ] dq- 1 



6[C;g]_ 







(15) 



These equations are second order linear differential equation with polynomial coefficients. The equation of Q[A; g]+ is 
inhomogeneous while the equation of 0[C;g] + is homogeneous. Our series solutions are particular solutions of these 
equations. The corresponding homogeneous equations are 



-(l-ff- 1 )^ 



dq 



dq 



Q[A;q]+ = 0, 



i s d 



dq 



dq 



6[C;g] + =0 



(16) 



These equations have one singular point q = of first kind and one q = oo of second kind. Our solutions are regular 
at q = 0. Because the singular point q = is first kind, we can consider their indicial equations. 



i/(v-l)+v = Q 



n(fi - 1) + fx = 



(17) 



This means that general solutions of these equations might have log divergence at q = 0. Polynomial approximation 
is not so bad because the convergent radii are infinite. However if we consider their logarithm, sometimes we are led 
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to wrong direction. In order to avoid it, we use differential equation though its validity in q >> 1 region is not so 
clear. 

Connected generating function W is log of grand partition function 9; W = —In 9 in Euclidean model. The 
expectation value of the matrix size N are defined as the derivative of W[0;/3] := — ln9[0;q] where the symbol O 
represents the type of sequence O = A,C; ip = (N)o = (d/d(3)W[0; (3\. The effective action r[0;<^] is defined by 
the Legendre transformation; T[0; ip] := W[0; /3] — f3<p. We will show their graph later. In the evaluation of matrix 




^ 10' 




FIG. 2: Grand partition functions and the expectation values of the matrix size iV and 27V for A and C, respectively. 




200 400 600 800 1000 1200 1400 1600 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

(JV) (N) 

FIG. 3: Effective potentials r[A]±, rL4]± as functions of the expectation value of matrix size, (-ZV). 
integrals, we omit the discussion on sign. Let us consider grand partition functions with sign; 

oo oo 

e[A;qU:=l+J22[AN-i](-q) N Q[C; q}- := 1 + £ Z[B N ](-)»- \ 2N ■ (18) 

N=2 N=l 
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We use the subscript ± in order to explain this difference. They follow to the following differential equations; 



d 2 d 
df +il + q ^ 



Q[A;q] 



d 2 . _ x . d 
d^ + iq + q ^ 



&{C;q}- = 2 



(19) 



Their graphs are included to the previous figure and we show their effective potential. The graphs for 0_ have cusp 
singularities. These effective potentials which are functions of the expectation value of the size of matrix are not single 
valued. There are two branches in SU and USp cases. These results imply the phase transition of these models with 
respect to the variation of the expectation value of the size of matrix. 



IV. CONCLUSION 



We reviewed non-Abelian Berry's phases in the USp matrix model as an effect of fcrmionic integration against the 
classical background. We encountered generalized monopoles. As such a monopole, we studied about the Tchrakian's 
five dimensional monopole. The charges of the Tchrakian's monopole and the Yang monopole are instanton number. 
The effect of Berry's phase is considered in the T-dualized model. In order to interpret them in terms of the original 
model, we performed matrix integrals. We obtain full result for all of classical gauge groups and their grand partition 
functions are computed. The grand partition functions are written as an exponential integral function and confluent 
hypcrgcometric functions. We showed graphs of grand partition functions of A and C series. In other words, matrix 
integral for the little IIB and little USp matrix model were exhibited. In these models, the expectation values of 
numbers N and their effective potentials were considered. There exist a problem of choice of sign. In the case of 
minus sign, we obtain cusp singularities. It shows that there are another phases of these model. From the figure the 
start points which correspond with q = are not stable. If we assume that these models have sense only when the 
matrix sizes are large enough, the reliable regions are around q — 1 and they are neighbor of those cusps. The graph 
for SU(N) show that the potential fall into the expectation value < N >= 1, while for USp(2N) the stable point 
looks like < N > = 0. We would like to compute matrix integral for full Itoyama-Tokura model. There are several 
interesting papers which might be related to our work -30j|. 
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APPENDIX A: BASIC TABLES ON CLASSICAL LIE GROUPS 



Let us list linear groups. We omit the Lorentz groups. 



{g€ M(N,K)\detg^0} , 



GL{N,K) 

U(N) := {g G M(N, 
SO(N) := {g G M(N, 

Jin 



1.9 1 




9 l = 9- 1 }, 



SL(N, K) := {g G M{N, K) | det g = 1} , 

SU{N) := U(N)DSL(N,C) , 
USp(2N) := {g G U(2N)\g t J 2N g = J 2N } , 



(Al) 



where K = M or C and N > 1. Group SO(l) is a discrete group which is isomorphic to {±1}- Group SU(1) is a 
unit group {1}. Groups U(l), SO(2) are isomorphic to each other. They are commutative and their fundamental 
groups are 7Ti([/(l)) ~ iri(SO(2)) ~ Z. We will omit these groups. Group SO(N)(N > 3) is not simply connected, 
but doubly connected; iri(SO(N)) ~ 7L 2 (N > 3). The universal covering group of SO(N) is denoted by Spin(N). 
Groups SU(N)(N > 2) are simply connected, group Spin(N) is formed in terms of Clifford algebra. Let us list their 
centers. We denote the center of a group G Z(G). The Schur's lemma show that the center of the general linear 
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group GL(N, C) is the group which consists of scalar operations Aljy(A G C x ) where k x = k \ {0}. 

Z(GL(N,K)) ~ A x , Z(SL(N,K)) ~ {z G K x ,z JV = 1} , 

Z(U(N)) ~ S 1 ~ {z e C x ||z| = 1} , Z{SU{N)) ~Z N ~{ze C x \z N = 1} , 

Z(SO(2N)) ~ Z 2 , Z(SO(2N + 1)) ~ 1 , 

Z{Spin{AN)) ~ Z 2 x Z 2 , Z{Spin{AN + 2)) ~ Z 4 , 

Z(Spin(2N + 1)) ~ Z 2 , Z(USp(2N)) ~ Z 2 , (A2) 

There are exact sequences; 1 — > 7Ti(50(iV)) — * Z(Spin(N)) — > Z(SO(N)) — » 1 (exact). Their Lie algebras are 
gl(N,K) := M (N, K) , sl(N,K) := {X G M(A, 2T)|TrA = 0} , 

u(iV) := {X G M(JV, C)| A* + X = 0} , su(A) := u(iV) n sZ(A, C) , 

so(N) := {X G M(JV,R)|X' + X = 0} , usp(2N) := {X G u(2A)|A* J 2A r + J 2A rX = 0} , (A3) 

where gl(N,C), sl(N,C) are complex Lie algebra and the remaining are real Lie algebra. The complexification of 
su(N) is sl(N,<C). We denote the complexification of Lie algebra so(N),usp(2N) so(N) c , usp(2N) c , respectively. 
Complex Lie algebra sZ(A,C)(A > 2),so(A) c (A > 3, A ^ 4), usp(2A) c (A > 1) are simple. There are several 
isomorphisms of Lie algebra; so(3) ~ su(2) ~ usp(2), so(5) ~ wsp(4) and so(4) ~ su(6). 

Let us focus on Lie algebra sl(N,C)(N > 2),so(N) c (N > 3, A ^ 4), usp(2N) c (N > 1). For shortening our 
notation, we use symbol sIn, son-, usp2N for them respectively. Let us consider their Cartan decomposition. For each 
algebra g, fj(g) denotes a Cartan subalgebra of g. 



f)(sZjv) = {diag(0i,02, • • • , fav), fa H h 0jv = 0} , t)(usp 2 N) = {[ . , A = diag(</>i, 2 , • 






Of course, they are not unique. t){usp2N) depends on the choice of J 2 at- A Cartan subalgebra () of a complex 
Lie algebra g is a maximal commutative subalgebra. The algebra f) acts on the algebra g by commutator; for 
each € f) and for each X G ad : f) — ► gZ(fl), ad(ff)(A) = [if, A]. The space g is a complex vector space. 
Every irreducible representation of commutative algebra is one-dimensional. Let V be such a irreducible part of 
g. The dimension of V is one and the action of ad(-ff) reduced to the multiplication of number cx(H) where the 
operation ad(H) is linear in H and these eigenvalues define a linear form on t); a £ \)* . Let us define a space 
V a := {X e g\[H, X] = a{H)X , for each H G f)} for each a G I)*. Let us define the root system with respect to 
this Cartan decomposition; $ := {a G dime /0,a^0}. Therefore g = (j Xae* ^a- Because g is a finite 
dimensional vector space, $ is a finite set. Let us list the root systems for sijv, so2N+i,usp2N, S02N- 

1. <&(sIn) = An-i case 

Sometimes the notation Ajy is used for the Lie algebra. We use this notation for the root sys- 
tem. Let us use the Cartan subalgebra of sijv which was introduced in the previous part; fj(sZjv) = 
{diag(0i, • • • , 4>n), fa + ' ' ' + <Pn = 0}. Let us define linear forms a^(i < j) such that aij(H) = fa — fa; 
(i < j,H G Fj(sZjv')). They and their counterparts form the root system $>(sIn) = Azv-i = 1 < * < j < 

A}. 

2. $(so 2JV ) = L>at 

Let us use the Cartan subalgebra f)(so 2 jv") = {0i J fa J • • • favJ}, where J = er 2 = [ 1 ] . Let us define 

linear forms a>ij, /3ij(i < j) such that aij(H) = <f>i — 4>j and (3ij(H) = <f>i + (f>j (i < j, H G ^(so 2 w+i)). 
The root system becomes $(so 2 at) = Dn = {±o%", iAj 1<*<J<A}. 

3. $(so 2A r +1 ) = Bjv 

Let us use the Cartan subalgebra f)(s02jv+i) = {<^i</ © 02^ © ■ ■ ■ 4>nJ ® (0)}- Let us define linear forms 
Oij,Pij,7k such that otij(H) = fa - fa, 0ij(H) = fa + <f> h 7fe (i?) = fa (1 < i < j < N, 1 < k < N, H G 

f)(s0 2 AT + i)). 

The root system becomes §{so2N+i) = Bw = {±a,j, ±/%, ±7fc 1<*<J<A, l<fc<A}. 
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4. <f>(usp2 N ) = C N 

Let us use the Cartan subalgebra f)(us£>2iv) = \ { * | i A = diag(^i, • • • ,4>n) \ ■ Let us define linear 

I \- iA / J 

forms aij,Pij^ k such that = fa- ^, = </>* + fa,lk{H) = fa (1 < i < j < N, 1 < fc < N, H G 

f)(so2JV+i))- The root system becomes $(usp2n) = CV = {ifly, ±/3y, ±27fe 1 < i < j < N, 1 < k < N }. 



1. Fundamental root system 



The following sets consist of positive roots with suitable ordering; 
{ctij, l<i<j<N} . 



1 A + 



2. B+ = {oij^ij,^, l<i<j<N,l<k<N} 



3. C 



N 



4. D 



N 



{otij, fa, 2 7fc 1 < i < j < N, 1 < k < N } . 
{ati3,0ij, l<i<j<N} . 



If we restrict the range of fa to real number, these roots define real linear forms. Therefore a root is denoted as a 
vector. For example «i2 of Dn can be expressed as (+1, —1,0, • • ■ ,0). However in An series, there are constraint 
fa + fa + ■ ■ ■ + fay = . This implies that any linear form have ambiguity caused byK = (1, 1, ■ ■ ■ ,1). Let us consider a 
vector (fa, ■ • • , fay) which satisfies the constraint. Let a be a linear form which maps such a vector into real number. 
Then a + Xk gives the same value as a, where A is an arbitrary real number. We would like to use the representation 
whose component satisfy the relation a\ + • • • + aj^ = 0. This constraint determines the ambiguous factor A. It might 
be a kind of a gauge slice. 

Now we can define the lexicographic order, i.e. for each real vector x = (x\,--- ,Xn), x > if and only if 
Xi = 0, • ■ ■ , Xk—i = and Xk > for some integer k. For two vectors x, y, x > y if and only if x — y > 0. This 
order is linear order of M. N such that every two elements can be comparable. In this sense our positive roots a satisfy 
the relation a > 0. Let $ be one of An, Bn,Cn, Dn and $ + be the positive root system. Because the ordering is 
linear, there is the minimal element s\. Again the set < I >+ \ Rsi has the minimal element. Let us find such a vector 
recursively, Sk+i = min{$ + \ Ks^} and we obtain a set {si, ■ ■ ■ , sn}- It is a basis of M. N . 

1. A N - s N = ai2, sjv-i = "23, • • ■ ,si = atf,N+i- 

2. B N : s N = ai2,sjv-i = "23, 

3. Cn- Sn = «12, SN-l = "23, 

4. Dn'- Sn = "12, Sn-1 = «23, ' ' ' , s 3 = Q!JV-2,JV-1) s 2 = Pn~1,N, Si = (XN-1,N 



,S2 

,s 2 ■ 

,S3 



ajv-i,iv, si = In- 

OtN-l,N> Si = 27jv- 
Q!JV-2JV-lj S2 = (3n-1,N,Si 



In fact, they are linearly independent. They form a lattice L(&) = ®^ =l 1si. At least in our case, L(Q) includes $. 
It implies that every root is a sum X^i a i s i- In these cases, the coefficients a% are positive integers for each positive 



root. For two vectors x = (x\, ■■ • ,xn),V = (yi, 
Let us define a matrix Cij = 2(sj, Sj)/(sj, Sj). 



,Dn) we obtain an inner product (x,y) = X\y\ 



xnVn- 



C[A 3 




C[B 3 




C[C 3 




C[D 3 



(2 


V 



o 

2 

1 -1 



-1 



(A6) 



These matrices are equal to Cartan matrices(We will omit the definition of Cartan matrix). 

Their determinants are det C[Ajv] = N + l,det C[Bn] = detC[CAr] = 2,detC[.Djv] = 4 (You can check these values 
with physical induction). They are orders of centers of corresponding simply connected groups. 



• A N : Si (1 < i < N), sj + Sj^! H h s l (1 < i < j < N) ; 

• B N : s t (1 < i < N), Sj + sj_i + • • • + s { (1 < i < j < N) , 
s 3 -\ h s l+1 + 2si H h 2si (1 < i < j < N) ; 
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• C N : s l (1 < i < N), Sj + Sj-! H h s 4 (1 < i < j < N) , 

2si H h 2s 2 +si (2 < i < iV), 

Sj H h s i+ i + 2si + • ■ ■ + 2s 2 + si (2 < i < j < N) ; 

• D N : Si{l<i< N), S] + H h s; (1 < i < j < N - 1) , 

Sj- + + • • • + s 3 + si (3 < j < N), Sj + Sj -_i + • • • + s 3 + s 2 + si (3 < j < A), 
Sj + Sj-! H h s i+x + 2sj H h 2s 3 + s 2 + si (3 < i < j < N) 

Weyl groups for A,B,C,D series are 

W[A N -i] = e N , W[B N ] = 6 N yi (Z/2Z) N , W[C N ] = & N x (Z/2Z) N , W[D N ] =S N yi (Z/2Z) N ~ 1 , (A7) 
APPENDIX B: VAN DER MONDE DETERMINANT 

In this section, we consider the measure on the Lie algebra consisting of Hcrmitian matrices. We start from the 
basics of Lie algebra sl(n). 

fl :={XeM(ra,C)|Tr(X)=0} , t) := {H e q\H = diag(/i 1( h 2 , ■ ■ ■ , h n )} . (Bl) 

Elements belonging to \) commute with each other. We define adjoint action of X £ g on Y G g: ad(X)(Y) := [X, Y]. 
For X, Y £ g, [ad(X), ad(K)] = ad([A, Y]). This is the adjoint representation of g. Especially, for iJi,iJ 2 £ h, 
[ad(-ffi), ad(i7 2 )] = and there are simultaneous eigenvectors T £ g of t), where T is not an clement in fj. Let T be 
such a generator, i.e. for arbitrary H £ [), there exists a complex number A such that [H, T] = XT. Suppose that for 
H 1: H 2 e f), the eigenvalues are Ai, A 2 £ C: [H U T] = AiT, [H 2 , T] = A 2 T. 

[oiffi + a 2 ff 2 ,T] = (aiAi +a 2 A 2 )T , (B2) 

This means that there is a linear map ^ from f) to the complex number field C: [H, T] — ip(H)T, H £ f). The dual 
space of f) is defined as f) v := Hom([),C). The linear map ip belongs to f) v . We call ip root. Let us define the set 
consisting of eigenvectors: s := {T £ g \ t)\3a £ f) v ,Vi? £ f), [i?, T] = a(H)T}. We define a set consisting of all roots: 
$ := {a £ f) v |[-ff, T] = a(H)T,T £ s}. Because g has finite dimensions, <£> should be a finite set. $ is called root 
system. Let a be a root. The eigenspace with respect to a is 

V a := {X £ q\VH £ (5, X] = a(H)X} . (B3) 

The dimension of this space is one. We define an inner product on g, k(X, Y) := Tr(ad(A)ad(F)). This inner product 
induces an inner product (■, ■) on f). 

(H, H) := Tr(ad(i?)ad( J ff)) = ^ a(H) 2 (B4) 

Because g is semi-simple Lie algebra, this inner product (H, H) is nondegenerate. We can construct orthonormal 
basis of this space: {e» }<=!,... >n , (ej,ej) = This inner product induces an isomorphism between f) and f) v . 

^:l)-[) V , <p(Hi)(H 2 ) :=(#!, ff 2 ) . (B5) 

Suppose that X is a Hermitian matrix. A measure on this space is defined as follows 

n n 

rj rj dSRXydsxy cx rj ^ rj dx Q (B6) 

i— 1 i<j z— 1 

A unitary matrix Z7 diagonalize the matrix X like X = [/AC/ -1 , where A is a diagonal matrix. Let us consider its 
differentials 

dX = diUAU- 1 ) = U{dK)U^ + (dU)AU- 1 + C/A(dL/ -1 ) , 

= U(dA- ad(A)(^)) U- 1 (B7) 

where the differential form A = U^ 1 (dU) takes value in the Lie algebra and the components along to the Cartan 
algebra do not affect. Therefore this form is defined on U(n)/T n , where T n is a maximal torus of U(n). The adjoint 
action of ad(A) yields the van-der Monde determinant term as the Jacobian: det | roo t(ad(A)) = Ila-root^' We 
can easily generalize this result to other groups. 
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APPENDIX C: EVALUATION OF INTEGRALS AND SPECIAL FUNCTIONS 



Here we present our evaluation of the residue calculus. 

A N _ 1 :fa = (N + l)/2-i, B N :<j>i=N-i + l, C N : fa = N - i + 1/2, D N : fa = N - i, (CI) 



where 1 < i < N. 
fa - <Aj = (j - i)- 



( N j-2 . . \ ( N j-1 

= h (n n j^j In n ^ I m 



! ^l = ll.2...( J -2 ) y ^±± 



AT! I H v ' I H J AT- AH 

\J=3 / \j=2 

fa = (.] -i),fa+ fa = 2N-(i + j). 

?m i 1 faff 2jV ~(W) ^ fa ^ ifaff 2^v- (, + ,,) 

/ n^^Vn^^i coo 



2 N ~ 1 N-N\ \ 11 2N-2j M 11 27V- ? 
\i=2 J / \i=2 

1 (2JV-3)!! 



2 N N-N\ (2AT-4)!! ' 
i>j = (j - i), fa + fa = 2(N + 1) - (i + j), fa = N — i + 1. 



, /i + 1) - (i + i) I l Jl 11 (27V + 3) - (i + i) 



1 (JV + l)-A [SL(N + l)-i 



. 1=1 



1 / -/t 2iV - j + 1 \ / JX- 2iV - 2j + 3 \ AT 



2^- TV! (2AT-2)!! 2N AT + 1 

1 (2AT-1)!! 
2 N + 1 N -M (2N-2)\\ 



(C7) 



zw ! _i (ffff (2jV + 2)-(z + J ) \ /**-3 (2JV + 2)-(i+j) , 



! I j2 2A^ - 2j + 2 J [n 2N - j + 2 j N + 1 

1 (2AT-1)!! AT + 1 AT 



(Cll) 
(C12) 
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<t>3 = U - *)> <j>i + h = (™ + !)-(* + j), 2cj> t = (27V + 1) - 2i 



Z[C 



N 



2 N N ■ TV! 



N j-1 

nn 

> 3=2 j=l 



(2N+l)-(i+j ) 
2N-(i+j) 



N j-1 



nn (2JV +!)-(» + j) 
llll (2iV + 2)-(i + i) 



n 



(27V + 1) - 2i 
2N~2i 



' AT 

n 



(27V + 1) - 2i 



l (27V + 2) - 2i 

Ifr-jj (2N + l)-(i + j) \ (y n 1 (2^1)-(i + j) 
2^7V • TV! I 1 = 1 1 = 1 27V - (i + j) J I 1 = 1 11 (27V + 2) - (< + j) 



'N-l 

n 



(27V +1) - 2« 



2N -2i 

N 



> N 

n 



(27V + 1) - 2i 
(27V + 2)-2i 



n 



27V - j 



2 N N-Nl \ 11 27V-2j + l 



JV 

n 



27V - 2j + 2 \ (27V - 1)!! (2N - 1)!! 
2iV - j + 1 J (27V -2)!! (27V)!! 



1 



(27V -2)!! TV (27V - 1)!! (27V - 1)!! 



2^ • TV! (27V - 3)!! 27V - 1 (27V - 2)!! (2TV)!! 
1 (27V -1)!! 



2 N + 1 N ■ N\ (27V -2)!! 

These numbers lead us to several special functions. The series 

sz_ s(s + 1) z 2 s(s + l)(s + 2) z 3 
<S>{s,t;z)-l + -- + t(t + 1) 2I + t(i + 1 )( t + 2 ) 3T + '" 

is called a confluent hypergeometric function. The product s(s + 1) ■ ■ • (s + n — 1) is denoted (s) n . The 
$(s, i; z) is a solution of the differential equation 



d 2 F 

' dz 2 



+ (7 - z) 



dF 

dz 



aF = 



(C13) 
(C14) 
(C15) 
(C16) 

(C17) 

(C18) 
(C19) 

(C20) 
function 

(C21) 



Ei(x) 



where C is the Euler constant; 



lim 



—dt- 



-dt 



= PV / — di , (.t > 0) 



C + In x + ^ 



k=l 



k-k\ ' 



C= lim | V i - Inn )= 0.577 215 664 90' 

re— >oo \ ^ — ' k j 
\k=l J 



The asymptotic behavior of the function Ei(x) with respect to the limit x — > +oo is 

\R n (x)\=0(\z\-»- 1 ) 



Ei(x) = — 
x 



,fc=0 



The series 



p-Fg(si) S2, ■ • ■ ,s p ] ti,t 2 , ■ ■ ■ ,t q ;z) — ^ 



fc=0 



(gl)fc(s2)fc • ■ • (Sp)k Z k 

(ti)fc(ia)fc- ••(<,)* fc! 



(C22) 
(C23) 

(C24) 

(C25) 
(C26) 



is called a generalized hypergeometric series. The confluent hypergeometric function is one of the generalized hyper- 
geometric series; ^(s, t; z) = i-Pi(s; i; z). There are several recursion relations and ladder operators. 



-$(a + 1,7+ l;z) = $(a + 1,7; z) - $(a,7;z) 

7 



-£-*(a, 7; = +1.7 +!;*)> (C27) 

dz 7 
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zv + a|$(a,7;z) = a$(a+l,7;z) (C28) 
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